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ABSTRACT 

It has been proposed that the cool cores of galaxy clusters are stably heated by cosmic 
rays (CRs). If this is the case, radio mini-halos, which are often found in the central 
regions of cool core clusters, may be attributed to the synchrotron emission from the 
CRs. Based on this idea, we investigate the radial profiles of the mini-halos. First, us- 
ing numerical simulations, we confirm that it is appropriate to assume that radiative 
cooling of the intracluster medium (ICM) is balanced with the heating by CR stream- 
ing. In these simulations, we assume that the streaming velocity of the CRs is the 
sound velocity of the ICM, and indicate that the heating is even more stable than the 
case where the streaming velocity is the Alfven velocity. Then, actually assuming the 
balance between cooling and heating, we estimate the radial profiles of CR pressure 
in six clusters only from X-ray observations. Since the CR protons interact with the 
ICM protons, we can predict the radial profiles of the resultant synchrotron radiation. 
We compare the predictions with the observed radial profiles of the mini-halos in the 
six clusters and find that they are consistent if the momentum spectra of the CRs are 
steep. These results may indicate that the cores are actually being heated by the CRs. 
We also predict broad-band spectra of the six clusters, and show that the non-thermal 
fluxes from the clusters are small in hard X-ray and gamma-ray bands. 

Key words: cosmic rays — galaxies: clusters: general — cooling flows — radia- 
tion mechanisms: nonthermal — galaxies: clusters: individual: A1835, A2029, A2390, 
Perseus, RXJ1347.5-1145, Ophiuchus 



1 INTRODUCTION 

Radio halos are diffuse radio emission that permeates clus- 
ters of galaxies l|Feretti fc Giovanninil 12008) . They are be- 
lieved to be synchrotron radiation from cosmic rays (CRs) 
in the intracluster medium (ICM). Although their surface 
brightness is low, they are often extended on scales of 
~Mpc. Since they are generally found in merging clus- 
ters, some phenomena associated with cluster mergers, such 
as shocks and turbulence in the ICM, may be responsible 
for the acceleration of th e CR particles fe.g. | Jaffd 



Roettiger. Stone, fc Burns! 19991: Takizawa fc Naitc 



Fuiita fc Sarazinl l200ll: lOhno. Takizawa. fc Shibata 



Fuiita. Takizawa. & Sarazin 2003; Brunetti et alj 
Brunetti fc Lazarianll2011 ). 




Among the radio-halos, so-called mini-halos are excep- 
tional. While they are small in number, they are found in 
the cool cores of cool core clusters, which have not been dis- 
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turbed by cluster mergers for a long time (e.g. lGovoni et ail 
120091 ; iMurgia et ail 120091 ). Their size (~ 100 kpc) is small 
compered with ordinary radio halos. The relation between 
the radio luminosity and the size is less clear for the 
mini- halos compared with the ordinary halos (|Murgia et al.l 
2009), whic h has made it difficult t o find the origin of the 
mini-halos. iGitti. Brunetti. fc Settil l|2002l ) proposed that 
the radio synchrotron emission comes from CR electrons 
reacce lerated by turbulence in the ICM. IPfromrner fc Enfilir] 
(|2004h argue that the electrons illuminating the mini-halos 
are of secondary origin; they are created throug h the inter- 
action between CR protons and ICM protons. iFuiita et al.l 
(2007) also proposed that the CR protons producing the 
secondary electrons are accelerated at shocks around the 
central active galactic nuclei (AGNs) in clusters. 

For cool cores, there is another mystery with the ICM. 
The radiative cooling time of the ICM in the cores is 
generally less than the age of the clusters, while that in 
the surrounding region is larger than the age. If there 
are no heating sources, the ICM in the core should cool 
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and cannot sustain the pressure from the surrounding. In 
this case, a gas flow towar d the c l uster centre should de- 
velop (a cooling flow; see Fabianl Il994l ). However, X-ray 
observations have shown that the mas sive cooling flows 
are not developing in cluster cores (e . g . Ikebe et al.l 19971; 
Makishima et al1l200ll; iPeterson et al.N200ll; iTamura et al.l 
200ll : iKaastra et all l200ll ; iMatsushita et al l 120021 ). ThusT 
some unknown heating sources prevent cooling flows from 
developing in the cores. 

In cool cores, acti ve AGNs interacting with the ICM 
are o fte n observed (e.g . iFabian et~a l. 2000; McNamara et all 



20001; iBlanton et al 
Mazzotta et all 2002 



20011; 



McNamara ct al. 



2002; Johnstone ct al 



. .Fujita et al 

2002; Kcmpner, Sarazin . fc Rickei 120021 ; iTakizawa et a 
200a : iFuiitaet al.l 120041 ). Through the AGN activities 



200 it 



CRs may be accelerated around the AGNs. The CRs may 
transfer the energy generated by the AGNs to the sur- 
rounding ICM ( e. g. [Tucker fc Rosnerlll983l; lRephaeh1ll987l; 



Pfrommer et"al] 120071 ; lJubelgas etahl 120081 ) 



20041 ; 

One of the tra nsfer mechanisms is CR st ream- 
ing dRephaelil Il979l; iBoehringer fc Morfilf I198 " 1 



Loewenstein, Zweibel. fc Begelmar] Il99ll ; iGuo fc OH 
2008), in which CRs interacting with Alfven waves move 
outwards in the cluster. The ICM is heated by the PdV 
work of the C Rs. We have shown th at the CRs can actually 
heat the ICM (|Fuiita fc Ohirall201ll . hereafter Paper I). The 
heating is fairly stable because it is not localised around 
the AGN, and the density dependence of the heating term 
is globally similar to the radiative cooling term. Moreover, 
the CR distribution is not much affected by the temporal 
change of the ICM distribution. 

If CR protons are heating cores, they should produce 
non-thermal emissions through their interaction with the 
thermal protons in the ICM. In that case, the radio mini- 
halos could be closely related to the heating of the cores. 
We calculated non-thermal spectra of cool core s that are 
heated by the CR protons (|Fuiita fc O hira 2012, hereafter 
Paper II). We compared the results with the radio observa- 
tions of the mini-halo in the Perseus cluster. We found that 
the momentum spectrum of the CRs must be steep so as to 
be consistent with the observations of the mini-halo. We also 
showed that the detection of non-thermal emissions from the 
CRs would be difficult in the hard X-ray and gamma-ray 
bands, because of the steep spectrum. Thus, it is desirable 
to find a way of studying the CR heating in more detail only 
in the radio band. 

In this paper, we show that CR heating can be stud- 
ied closely by the radial profiles of mini-halos. If the CRs 
are stably heating a cool core, radiative cooling must be 
balanced with the heating at each radius. Assuming they 
balance, we can determine the radial profile of the CR pres- 
sure using only X-ray data for a cluster. Once we obtain 
the CR pressure profile, we can predict the radial profile of 
synchrotron radiation from the CRs and we can compare 
the predicted profile with the observed one. The idea of us- 
ing the balance between the radiative cooling and the CR 
heatin g has been proposed bv lColafrancesco fc Marchegianil 
(2008). They considered Coulomb and hadronic interactions 
as the CR heating. However, since it has been indicated 
that CR streaming is a much more effective heating source 



|Guo fc Ohll2008l . Paper I), we consider the CR streaming 
as the heating source in cores. 

This paper is organised as follows. In Section [5] we de- 
scribe our model on the estimation of CR pressure and the 
non-thermal emissions from the CRs. In Section[3] we check 
whether the CR heating is actually balanced with the radia- 
tive cooling using the results of numerical simulations. In 
that section, we also show that the CR heating is even more 
stable if the streaming velocity of the CRs is the sound ve- 
locity of the ICM, compared to the case where the streaming 
velocity is the Alfven velocity. In Section [4] we predict radial 
profiles of mini-halos and compare them with the observa- 
tions. Finally, Section[5]is devoted to discussion and conclu- 
sions. Throughout this paper we assume a ACDM cosmology 
with H = 71 kms" 1 Mpc -1 , f2 ro = 0.27, and Q A = 0.73. 
We consider protons as CRs unless otherwise mentioned. 



Rephaeh fc silkT 1 19951; IColafrancesco. Par, fc De Ruiulal 2 MODELS 



2.1 CR Pressure 

If the heating via CR streaming is balanced with radiative 
cooling in a cool core, it must be 



H st » n e A(T, Z) . 



(1) 



where n e is the electron density of the ICM. The heating 
rate of the CR streaming is given by 

dP c 



H, 



~V B f 



dr 



(2) 



where v st is the streaming velocity of the CRs, and P c is 
the CR pressure (Paper I). In our models, dP c /dr < 0. We 
approximate the cooling function A by 



A(T, Z) = 2.41 x 10" 



+1.39 x 10" 



ii> + U.I [ — 



0.02 + 0.1 



(!) 



erg cm 



(3) 



where T is the temperature and Z is the metal abundance 
of the ICM. This function a pproximates the one derived by 
ISutherland fc Dopital (|l993l ) for T > 10 5 K and Z < 1 Z Q . 
From equation (O, we obtain CR pressure: 



Po(r) 



njA 

Vst 



dr . 



(4) 



where r max is the radius where P c approaches zero. Since nj! 
rapidly decreases outside the core, the integration is insen- 
sitive to the value if we assume that it is much larger than 
the core radius. Thus, we assume that r max = 1 Mpc. 

CRs move outwards as a whole in the cluster with waves 
that scatter the CRs. If the waves are Alfven waves, the 
streaming velocity v B t may be the Alfven velocity va- How- 
ever, it has been indicated that the CR streaming velocity 
may be much larger than the Alfven velocity in hot ICM, 
because in hot plasma the Alfven waves may suffer strong 
resonant damping at small wave lengths by thermal pro- 
tons. In this case, the soun d velocity c s may be appropriate 
as the streaming velocity (|Holman. Ionson. fc Scott] Il979l ; 
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Enfili n et al.l l201ll ). In this paper, we mainly consider the 
case of v Bt = c s . 

In this study, we do not discuss the injection of CRs into 
the ICM in detail. They may be directly injected at shocks 
dFuiita et alj|2007h. or th ey may be transported by buoyant 
bubbles fc Ohll2008l ). In the latter case, the bubbles are 
created through AGN activities, which also accelerate CRs 
around the bubbles. The bubbles contain the CRs and carry 
them from the AGN out to large distances. As the bubbles 
adiabatically rise, the CRs may escape from the bubbles into 
the ICM or they may be injected into the ICM through the 
shredding of the bubbles by Rayleigh- Taylor and Kelvin- 
Helmholtz instabilities. While the adiabatic expansion re- 
duces the energies of CRs, it does not change the slope of 
the momentum spectrum. Observations have shown that the 
bubbles with a size of ~ 10 kpc are often found at the cen- 
tre of cores, while large scale (~ 100 kpc) shocks appear t o 
be rare (e.g. iBirzan et afll2004l : iMcNamara fc Nulsen|[2012T ) . 
This may suggest that CR injection via bubbles is more 
common. 



2.2 Synchrotron Emission from CRs 

Although we can obtain P c (r) from equation Q for given 
temperature and density distributions and the streaming ve- 
locity, we need to specify the momentum spectrum of the 
CRs to calculate the synchrotron spectrum of the CRs. For 
that purpose, we adopt the following spectrum: 



N(p, r) = A„(r)p x e 



-P/Pn 



(5) 



where p is the CR momentum, x is the index, and p max is 
the cutoff momentum. 

The hadronic lo ss of CR protons could modify the slope 
of the spectrum (e.g. lBrunetti et al . 2004; E nfilin et al, 2007; 
iGuo fc Ohll2008l ). In our case, the escape time of CRs from 
the core (~ 10 8 -10 9 yr) is smaller than the cooling time 
via hadronic loss. Thus, we ignore the modification of the 
spectrum by the loss. 

We assume that x and p ma x do not depend on r (Pa- 
per II). On the other hand, the normalisation of the spec- 
trum, A cr (r), depends on r and it is determined from the 
following equation for given P c (r): 



Mr) = 3 



P 2 N(p,r) 
p 2 + m 2 c- 



:dp . 



(6) 



We determine the cutoff momentum p max assuming that the 
high-energy CRs are accelerated at shocks formed through 
the activities of the central AGN (see Paper II). The cut- 
off momentum is determined by the age of the shock, and 
/-'max'-- 10 5 TeV if the energy input through an explosive 
activity of the AGN is ~ 10 60 erg. Although the actual accel- 
eration mechanism is not clear (e.g. CRs may be accele rated 
nearby the black hole and form jets; ISikora e t al. 2005J), the 
following results are not sensitive to p max - The minimum 
momentum, p m in, will be set in Section 13.11 On the other 
hand, as will be shown later, the synchrotron spectrum and 
luminosity are sensitive to the index x. 

Once we fix the momentum spectrum of the CRs at 
each radius, we can calculate synchrotron and other non- 
thermal emissions from secondary electrons that are gener- 
ated through the interaction between CR protons and ICM 



protons, and through the decay of charged pions. We can 
also calc ulate 7r°-decay g amma -rays. We use the radiation 
model of iFang fc Zhanel (|2008T> and the p roton-proton col- 
lision model of iKarlsson fc Kamaei (|2008l ); the details are 
described in iFuiita et all (|2009h . 



3 COMPARISON WITH NUMERICAL 
SIMULATIONS 

3.1 Balance between Heating and Cooling 

Before we study real clusters, we investigate whether the 
results of numerical simulations are well approximated by 
equation (TTj). We calculate P c {r) based on the models in 
Papers I and II, although we make a few modifications as 
follows. We adopt the gravitational potential and the ini- 
tial ICM profile of Model LCRs in Paper II. This model is 
originally constructed to reproduce the observations of the 
Perseus cluster. We do not include thermal conduction, and 
the CR diffusion can be ignored. The CR streaming velocity 
Vst is the sound velocity c s ■ We reevaluate the parameters for 
Coulomb and hadronic losses of CR protons. This is because 
although we found that the index must be large (x ~ 3) in 
Paper II, the model we used to obtain P c {r) in Paper II 
did not include the effect. Th e parameters of t he losses are 
calculated as follows (see also lGuo fc OhlfeoOSl ). 

The energy loss of a proton with velocity of v = /3c b y 
Coulomb loss is given bv lMannheim fc Schlickeiserl (| 19941) : 

P 2 



( dE v 
\ dt 



-4.96 x 10" 



V cm" 



■ erg s 



(7) 



/3 3 +a4 

where x m = 0.0286 (T/2 x 10 6 K) 1/2 . For the momentum 
spectrum (|5]l. the overall Coulomb loss rate is 



n(** 

V dt 



dp 



= -7.3 x 10" 



erg cm" 



-1 -3 
erg s cm 



(8) 



where e c is the CR energy density. In this subsection, we as- 
sume that x = 3 and T = 8 keV, which are the typical values 
for the clusters we later investigate. The minimum momen- 
tum is set to be p m i n c = 137 MeV, which corresp onds to the 
kineti c energy of 10 MeV and was adopted by IGuo fc Ohl 
(2008). Around this momentum, the Coulomb loss of a pro- 
ton is maximum (/3 ~ x m in equation [7]). 

The energy-loss rate of a CR proton due to pion pro- 
duction is approximately given by 

fdE p \ 
V dt )h 

where njv is the nucleon density, a vv is the pp cross-section, 
Pthr is the threshold momentum, T v is the k inetic energy 
of th e proton, and 9 is the step-function (|Enfilin et al.l 
| 2007l). For un, &vv , and pthn we adopted the ones used in 
lEnfilin et alj l|2007l ). The overall hadronic loss rate is 



-0.5n N a pp cTp6(p - p t hr) , 



(9) 



r„ = 



N 



( dE v 
\ dt 



dp 



-1.5 x 10" 
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and (b) density profiles for Model 1. 
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Figure 2. Profiles of the ratios P c /Pg (solid) and Pb jP a (dotted) 
at t = 9 Gyr for Model 1. 



erg cm - 



-1 -3 

erg s cm 



(10) 



Thus, the energy loss rate is ri oss = — r c — I\ = £ c n c e c , 
where ( c = 7.5 x 10" 16 cm 3 s -1 . Since ~ 1/6 of the in- 
elastic en ergy goes into secondary electrons during hadroni c 
collisions l|Mannheim fc Schlickeiseilll994l ; lGuo fe Ohf2 008). 
the heating rate by the Coulomb and hadronic collisions is 
H coil = -Fc-Fk/Q = r/ c n c e c , where r) c = 7.4x 10~ 16 cm 3 s _1 . 
We use these values in equations (3) and (4) in Paper II. 
The values are no t much different from the ones derived by 
iGuo fc Ohl | 20081 ) and adopted by us for x = 2.4 using a 
slightly different formula (£ c = 7.51 x 10~ 16 cm 3 s _1 and 
rj c = 2.63 x 10" 16 cm 3 s _1 ). 

In this subsection, we consider two models, which are 
based on Model LCRs in Paper II. Contrary to Model LCRs, 
we adopted the above £ c and r\ c - Moreover, we assume that 
the adiabatic index of the CRs is y c = 5/3 instead of 4/3 as 
an extreme case, because the index approaches 5/3 if most 




r (kpc) 

Figure 3. CR pressure profiles for Model 1. 




r (kpc) 



Figure 4. Relative importance of CR streaming (H B t) and colli- 
sional heating (H co u) at t = 9 Gyr for Model 1. 



of the CRs have low energies. In Papers I and II, the total 
CR injection rate by the AGN was given by I/agn = eMc 2 , 
where M is the inflow rate of the gas toward the cluster 
centre. If the CRs are gradually injected into the ICM as 
the buoyant bubbles rise, their injection should not be lo- 
calised. Thus, we give the CR injection rate per unit vol- 
ume by S c oc r~" for 20 < r < 150 kpc. The lower bound 
is ba s ed on the observ ations of the bubbles (|Birzan et al.l 
120041 ; IGuo fc Ohl 120081 ). Although the upper bound is un- 
certain, the results are not sensitive to the value, because 
r 3 S c is a decreasing function of r if v > 3. In one model, 
we assume that the efficiency is e — 2.5 x 10~ 4 , and the 
radial dependence is v = 3.1. These values are the same as 
those in Model LCRs in Paper II. We refer to this model as 
Model 1. We consider another model, which is the same as 
Model 1 but the AGN feedback is stronger (Model 2). We 
assume that e = 2.5 x 10~ 3 , and that the energy input is 
more centrally concentrated (y = 3.5). We have shown that 
the CR heating is locally unstable with these parameters (e 
and v) if the streaming velocity is the Alfven velocity and 
if thermal conduction is not included, because the heating 
is too strong and too centrally concentrated (Paper I). We 
assume that B = 10 (n e /0.016 cm" 3 ) 2/3 ^G, which is the 
same as that in Model LCRs in Paper II. 
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r (kpc) 



Figure 5. Same as Figure[T]but for Model 2. 




r (kpc) 

Figure 6. Same as Figure[2]but for Model 2. 



It is to be noted that M represents the mass flow that 
passes the inner boundary (r = 5 kpc), and that not all the 
mass is swallowed by the central black hole. Most of the gas 
may be consumed by star form ation in the central galaxy 
(e.g. IMcNamara fc Nulserj|20l3 ) . For example, if 2.5% of M 
goes to the neighbourhood of the black hole, and 10% of 
the rest mass energy is converted to CRs, the efficiency is 
e = 2.5 x 10" 3 . Although the conversion rate to CRs (10%) 
may be rather high, recent studies have shown that most 
of the AGNs at present se em to be in radiatively ineffi cient 
mode or radio mode (e.g. IMcNamara fc Nulsen I2012T ). and 
that the inertia of AGN jet s seems to be dominated by CR 
protons l|Sikora et alj|2005h . Thus, the conversion rate may 
not be unrealistic. The injected CRs eventually stream out 
of the cluster, to beyond the virial radius. 

Figures HHH are the results for Model 1. Figures [1] 
and [2] are almost identical to the corresponding figures for 
Model LCRs in Paper II (Figures 8 and 9 in that paper). 
The cluster is isothermal and in hydro-statistic equilibrium 




r (kpc) 

Figure 7. Same as Figure \3\ but for Model 2. 




r (kpc) 

Figure 8. Same as Figure [4] but for Model 2. 



at t = (Figure [TJ. The mass inflow rate gradually in- 
creases with time. At t = 12 Gyr, the mass inflow rate 
is M = 144 Mq yr . If there is no heating, it should 
be 760 Mq yr -1 (Paper I). In Figure [2j we show the ra- 
tios of the CR pressure P c and the magnetic pressure Pb 
to the ICM pressure P g . Figure [3] shows the evolution of 




200 



400 200 
r(kpc) 



400 



Figure 9. Profiles of synchrotron radiation when we use P c {r) di- 
rectly calculated via numerical simulations (solid lines) and when 
we use Pc{t) estimated by equation ||4} (dotted lines), (a) Model 1 
and (b) Model 2. The redshift of the cluster is set to be z = 0.0179 
(the Perseus cluster). 
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r (kpc) 



Figure 10. Same as Figure ITl but for Model 3. 




r (kpc) 

Figure 11. Same as Figure |21 but for Model 3. 



P c (r). For r > 100 kpc, the density of CRs nearly follows 
p c oc P}^ lc oc r~ 2 , or P c oc r -10 / 3 . This is because of the 
continuity, almost constant sound velocity (streaming veloc- 
ity), and inefficient cooling in the outer region of the clus- 
ter. However, the slope is less steep for r < 100 kpc, be- 
cause CRs are injected not at the cluster centre but mainly 
at 20 < r < 150 kpc. Figure [4] demonstrates that the CR 
heating, H Bt , and radiative cooling, n 2 A, are well balanced 
except for the innermost and the outermost regions. The fig- 
ure shows that 0.5 < # st /(n 2 A) < 0.8 for 20 < r < 200 kpc. 
Figures [SHE] demonstrate the results for Model 2. The CR 
heating is stable in spite of the large e and v. This means 
that the CR heating with v st = c 3 is even more stable than 
that with v s t = va- Because of the strong feedback, radia- 
tive cooling is almost cancelled out. At t — 12 Gyr, the mass 
inflow rate is only M = 17M0 yr _1 . The reason of the addi- 
tional stability is that if the ICM temperature decreases at 
the cluster centre, the streaming velocity v st = c s decreases 
there. Thus, the escape of the CRs from the centre delays, 




r (kpc) 

Figure 12. Same as Figure[3]but for Model 3. 




r (kpc) 

Figure 13. Same as Figure|4]but for Model 3. 

which increases \dP c /dr\ and the heating rate (equation [5]). 
Figure H shows that 0.7 < H st /(n 2 e A) < 1 for r < 100 kpc, 
and H st /(nlA) ~ 0.5 at r ~ 200 kpc. The well-established 
balance between heating and cooling both for Models 1 and 2 
suggests that the ICM evolves in such a way as to achieve 
the balance regardless of the details of the AGN heating 
model. In other words, equation (fl) is the most important 
factor that determines the ICM profile in the core. It is to be 
noted that the bulk of CR energy is dissipated as heat within 
the core. In Model 2 for example, the CR energy injection 
rate from the AGN is 2 x 10 45 ergs" 1 at t = 12 Gyr. At the 
same time, the flux of CR enthalpy escaping the core (say, at 
r = 300 kpc) is only 47T7 c u st r 2 P c /(7 c - 1) ~ 3 x 10 44 ergs -1 . 

Figure shows the radial profiles of synchrotron radio 
emission for Models 1 and 2 at t = 9 Gyr. We chose x = 
2.9, because the value is favourable when compared with 
observations (see Section T4.2p . We compare the profiles that 
are calculated using P c {r) in Figures [3] and [7] with those 
using P c (r) derived from n e (r), T(r), and v st (r)(= c s (r)) 
by equation Q . Figure [5] shows that their differences are 
small in the core region (r < 200 kpc) , which means that the 
approximation by equation {TJ is good. However, outside the 
core region, the results using equation Q is overestimated 
by a factor of a few, because the CR heating is not balanced 
with radiative cooling (Figures [4] and [8]). 
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Table 1. Cluster profiles 



Name 



Profiles' 1 



Rcfs 



Perseus 



A1835 



A2029 



A2390 



RXJ1347 



0.0179 n e (r) 



0.2532 



0.0765 



0.2280 



0.4510 



Ophiuchus 0.028 



l + (r/18) 3 il+^VSS) 5 ] 175 [l+tr/m) 2 ] - 87 

T(r) - 7 Q ^(r/ro) 3 

J " _ '• U 2.3+(r/70) : * 
Z(r) = 0.46 + 0.27 exp 

n e (r) 
n e (r) 



81, 



0.42 exp 



0.(11)2 



[l+(r/41)j^ ( fOT T < 231 ) 

[i+tXV-^ (for r > 231) 



T(r) = 7.8 - 3.7exp 
Z(r) = 0.28 + 0.15 exp 



.1 (JLY 

2 \4SJ 



2 \57J 
2 (r/85.5)- 1164 



n e (r) 2 - (0.0168) 2 [1 L r . /85 . 5) a ] i.os [i + ( r /923)^]" SBei 
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a The units for n e , T, Z, and r are cm 3 , keV, Zq, and kpc, respectively 
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Allen. Schmidt, fc Fabianl j2002f l 



Gitti, Piffaretti, fc Schindlerl ||2007[) 



lOta et al.l (|2008h 



J iNevalainen et al.l ||200I 
k iFuiita et ah! |2008) 



3.2 Structure of the magnetic field 



In our model, CRs are scattered by Alfven waves. We 
have implicitly assumed that the gradient of P c and 
the background magnetic field are parallel. However, 
magne tic field is o ften thought to be tangled in the 
ICM (Sarazin 1986), which may affect the CR stream- 
ing. It has been indicated that even if the gradient of 
P c and the background magnetic field are not paral- 
lel, t he plasma i s unstable and waves are excited by 
CRs (iBelll 120051; iRakowski. Laming, fc Ghavamianl 120081 : 



iBvkov. Osipov. fc Ellison! 1201 ll : iRogachevskii et al.l l2012h 
Thus, we expect that CRs are scattered by the waves. How- 
ever, the effective streaming velocity for a tangled magnetic 
field is uncertain, and previous studies have not dealt with 
this problem (e.g . lLoewenstein. Zweibel. fc Begelmanlll99ll ; 
iGuo fc Ohl 12003 ). Although the detailed study about it is 



beyond the scope of our study, we investigate the effect in a 
simple way. 



We assume that the streaming velocity for randomly 
tangled magnetic fields is reduced by a factor of a from the 
one for the radially stretching magnetic field (v st — ac a ). 
For a given magnetic field line, forward or backward Alfven 
waves can exist. However, when the waves are excited by 
CRs, they are driven only in the direction closer to the CR 
pressure grad ient , or outwards in a cluster in our case (e.g. 
lLongairl Il994f ). Owing to the scatter by the waves, a CR 
particle is not bound by a magnetic field line. On average, 
CRs move outwards along with the waves. If the angle be- 
tween the magnetic field line and the direction of the CR 
pressure gradient is 8 71 -/2), the effective wave velocity 
in the direction of the CR pressure gradient may be smaller 
than the sound velocity by a factor of cos9, and a may be 
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Figure 14. Ratio of CR pressure P c to ICM pressure P g . 
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Figures I10H13I show the results for the same parameters of 
Model 2 but for a = 0.5. We refer to this model as Model 3. 
At t = 12 Gyr, the mass inflow rate is M = 16 Mq yr _1 . 
The results are not much different from those of Model 2. 
However, the smaller streaming velocity makes the escape of 
CRs slower and increases P c (Figures [71 and [12)) . It appears 
that the increase of P c cancels out the decrease of v st in 
equation ((2|). Since the difference between Models 2 and 3 
is not significant, we assume that a = 1 in the following 
sections, unless otherwise mentioned. 



4 COMPARISON WITH OBSERVATIONS 
4.1 ICM models 

We study six clusters (Perseus, A1835, A2029, A2390, 
RXJ1347, and Ophiuchus) for which the rad ial profiles of 
radio mini-halos have been studied in detail bv lMurgia et al.l 
(2009, see also the references in the paper). We need to find 
n e (r), T(r), and Z(r) for the ICM of these clusters to obtain 
Pc(r) using equation Q. In Table 1, we show the ICM pro- 
files we adopted. Some references give the functional forms. 
For other references, we fit the observational data with ap- 
propriate functions. 



4.2 Results 

Using the ICM models in Section T4. II we can estimate P c (r) 
by equation Q for v s t = c s . Figure [TU shows the ratio of CR 
pressure P c to ICM pressure P g for the six clusters. The ICM 
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Figure 16. Profiles of synchrotron radiation of the Perseus clus- 
ter at v = 1.4 GHz. (a) Streaming velocity is Alfven velocity, vj±. 
The solid line shows x = 2.9, and the dotted line shows x = 3.5. 
(b) B = 3 (n e /0.016 cm" 3 ) 2 / 3 fj,G. The solid line shows x = 2.9, 
and t he dotted line shows x = 2.8. Filled cir cles are the observa- 
tions jPedlar et al.lll99fjl : iMurgia et alll2009h 



pressures are derived from the ICM models in Section 14.11 
The CR pressures at the cluster centres are less than 10 %. 
In Figure [15] we show the predicted surface brightness pro- 
files of synchrotron radiation for P c shown in Figure [T?] For 
these calculations, we need to specify B(r). In Figure [151 we 
assume that B = 10(n e /0. 016cm~ 3 ) 2 / 3 /jG. For comparison, 
we show the observations (IMurgia et al.ll2009l ; iPedlar et al.l 
1 19901 ; iBacchi et al1l2003l ; iGitti et al.ll2007l V As can be seen, 
the observations are consistent with x ~ 2.9. The observa- 
tions of all the six clusters are reproduced fairly well under 
the simple assumption of equation {1} and almost common 
index x ~ 3. In particular, the predicted slopes of the radial 
profiles are consistent with the observations, and the slopes 
are less uncertain than x. As was discussed in Section [3.11 
the surface brightness may be overestimated a factor of few 
outside the core (r > 200 kpc). Even so, the agreement be- 
tween the predictions and the observations is good. Among 
the six clusters, the consistency between the prediction and 
the observation s of A 2390 may be the worst (Figure [Tot O . 
IVikhlinin etafl (|2006h indicated that the central region of 
the cluster is strongly disturbed by the central AGN, which 
might affect the results. 

Figure [16] is the same as Figure 115b (Perseus) but for 
v s t = va (Figure 116b .) and for a reduced magnetic field 
strength (Figure 116b ). In general, the Alfven velocity is 
smaller than the sound velocity in the ICM. Thus, when 
v 3 t = Va, the escape of the CRs from the core delays and 
P c becomes larger. Accordingly, for a given spectral shape, 
the synchrotron emission becomes larger than that in the 
case of Vst = c s (the solid line in Figure 116b and that in 
Figure ITSb) . However, if we increase the spectral index to 
x — 3.5, the observations can be reproduced (the dotted line 
in Figure 116b ) . Although we do not present the result in the 
figure, we have also studied the profile for a reduced stream- 
ing velocity (v st = 0.5c s ). As is discussed in Section 13.21 
the reduced velocity increases P c . However, the effect is not 
significant, and the index of x = 3.0 can reproduce the ob- 
servations. In Figure 116b . we present the results when we 
reduced the mag netic field to B = 3 (n e /0.016cm _3 ) 2/3 fj,G. 
The synchrotron emission reduces only a factor of a few for a 
spectral index of x ~ 3 (the solid line in Figure [T&b and that 
in Figure [Tab ) . If we slightly reduce the index to x — 2.8, 



the observations are reproduced again (the dotted line in 
Figure [T&b ). Thus, the uncertainty of magnetic fields does 
not much affect the results. 

In Figure [T71 we present the broad-band spectra of the 
six clusters when v st = c s ■ We calculated non-thermal emis- 
sions from secondary electrons created through the decay of 
charged pions that are generated through the interaction be- 
tween CR protons and ICM protons. We also calculated 7r°- 
decay gamma-ray emissions. We did not include non-thermal 
emissions from primary electrons that are directly acceler- 
ated in the ICM, because we have no information on them. 
We adopted x = 2.9 because it can fairly reproduce the ra- 
dio profiles for the clusters (Figure [T5"1l . Figure [T7l shows that 
except for the Perseus cluster, GeV gamma-ray fluxes are 
much smaller than 10 -12 erg cm" 2 s" 1 , which is the typica l 
current upper limit set by Fermi ( Ac kermann et al"1l2010l ). 
For the Perseus cluster, the contamination of the central 
AGN increases the upp er-limit to ~ 10 -11 erg cm" 2 s -12 
l| Ackermann et al"1l2010l ). Because of the steep CR momen- 
tum spectra, the fluxes in the hard X-ray band and the 
TeV gamma-ray band are small. Thus, the detection in these 
bands would be very difficult, which i s consistent with re - 
cent TeV gamma-ray observations (e.g. Al eksic et al.ll2012l ). 
If non-thermal emissions are observed in these bands in some 
clusters in the future, they might come form primary elec- 
trons. 



5 DISCUSSION AND CONCLUSIONS 

We have shown that the radial profiles of radio mini-halos 
observed in galaxy clusters can be explained by the syn- 
chrotron emission from the CRs that heat the cool cores of 
the clusters. In our model, the cores are being heated via 
CR streaming. 

First, we investigated whether the CR heating is bal- 
anced with radiative cooling of the ICM using numerical 
simulations. We assumed that the streaming velocity of the 
CRs (v st ) is the sound velocity of the ICM (c s ), because it 
has been pointed out that it is more appropriate than Alfven 
velocity (va) in hot ICM. Moreover, we found that the heat- 
ing is more stable when v st = c s than when v st = va- We 
showed that the CR heating is actually balanced with the 
radiative cooling inside the core, although the heating falls 
short of the cooling a factor of a few outside the core. 

Assuming that the CR heating is balanced with the 
cooling in the cores of six real clusters (Perseus, A1835, 
A2029, A2390, RXJ1347, and Ophiuchus), we estimated the 
radial profiles of CR pressure only from X-ray observations. 
The results showed that the CR pressure at the cluster cen- 
tres are less than 10 %. Assuming that the momentum spec- 
tra of the CRs are given by a power-law and the CR stream- 
ing velocity is the sound velocity, we calculated radial pro- 
files of the synchrotron emission from secondary electrons 
created through proton-proton interaction. We found that 
the profiles match the observations if the index of the CR 
momentum spectra is x ~ 3. If the CR streaming veloc- 
ity is the Alfven velocity, the index must be x ~ 3.5 to 
be consistent with the observations. Uncertainties of mag- 
netic fields in the ICM do not much affect the results. We 
have compared our predictions with the observations of only 
six clusters, and our model has free parameters such as x. 
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Figure 17. Broad-band spectra of the six clusters. Synchrotron 
radiation (dotted line) , inverse Compton scattering off cosmic mi- 
crowave background (solid line), and non-thermal bremsstrahlung 
(dashed line) are of the secondary electrons. The 7r°-decay 
gamma-rays are shown by the two-dot-dashed line. For compar- 
ison, the thermal bremsstrahlung from the ICM is shown by the 
dot-dashed line. Rcdshifts have been corrected. 



However, if radial profiles of mini-halos are studied for many 
other clusters in the future, our model can be examined more 
extensively through the comparison with the observations. 

We also estimated broad-band spectra for the six clus- 
ters when x ~ 3. Owing to the steep spectra, hard X-ray 
and gamma-ray fluxes are small, which makes it difficult to 
detect the non-thermal emissions from the CRs in the clus- 
ters in those bands. In the radio band, future observations 
of the spectral index will give us information on the momen- 
tum spectra of CRs. For the Perseus cluster, we have shown 
that our predicted index is consistent with the observations 
(Paper II). 
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